Neutrino fluxes from CNO cycle in the Sun in the 
non stationary case with mixing. 
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Abstract. The computational analyses is presented of the non stationary case with 
mixing of the solar model when the neutrino flux F13 from the decay of 13 A is higher 
\ than a standard solar model predicts. 
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1. Introduction 

Recently it has been argued PQ- [3] that to achieve agreement of the solar surface heavy 
element abundance with the results of helioseismology [H [5] probably one needs to 
introduce some non stationary process in the solar model. In the paper [6] we addressed 
the question how the mixing in the Sun can change the fluxes of neutrinos generated in 
the CNO cycle. The point which drew our attention was that there is a large difference 
of abundance of carbon in the central and peripheral layers due to burning of carbon 
in thermonuclear reactions in the center of the Sun. This makes possible to obtain a 
relatively large increase of the flux F13 from the decay of 13 N even for a very mild mixing 
between the central and peripheral layers of the Sun while all other neutrino fluxes and 
other observable are not practically changed. Looking for the effect of mixing we have 
omitted coordinates in our consideration neglecting the specific character of the mixing, 
paying attention only on the evolutionary part of equations. 

' dX( 12 C)/dt = -X( 12 C)X( 12 C) + ^\{ 15 N)X( 15 N)+ 
+k(X ( 12 C)-X( 12 C)) 
dX( 13 C)/dt = -\( 13 C)X( 13 C) + if \{ 12 C)X{ 12 C) - kX{ 13 C) 
* dX( 14 N)/dt = -\( 1A N)X{ 14 N) + ^\{ 13 C)X{ 13 C)+ [) 

+k(X ( li N) - X^N)) 
^ dX( 15 N)/dt = -\{ 15 N)X( 15 N) + ^X( 14 N)X( 1A N) - kX( 15 N) 

where X( A Z) and Xq( a Z) - the abundance of isotope A Z in the central and peripheral 
regions, \( A Z) = pNaX i 1 H) < A Zp > - the (p + A Z) reaction nuclear rate [3 E], 
k - the mass transport coefficient. Here the mixing is understood as a spherically 
asymmetric process when the masses coming to the core of the Sun and leaving the 
core go by different paths. So this differs from gravitational settling and diffusion and 
in this sense is a 3D extension of the standard solar model which is known to be a ID 
model. But we followed a simplified approach. We used a weak mixing (which means 
that it involves very weak mass transport) which proceeds fast in comparison with the 
rates of thermonuclear reactions in the Sun. This enabled us to divide the internal zone 
(0.6 solar radius) of the Sun on spherical layers (onion like structure) isolated one from 
another. The set (CQ) has been solved for each spherical layer with the initial conditions 
of the standard solar model at the moment when the mixing has been initiated. Of 
course the results obtained in this way may serve only as guidance. Our aim was to 
get the averaged values and to see how they differ from the ones given by the standard 
solar model. Looking for the solutions of this set of equations we were particularly 
cautious not to go beyond the limits for the parameters of the solar model fixed by 
helioseismology. This remarkable selectivity of the mixing to the neutrino fluxes (_F 13 
is increased while F15 from the decay of 15 is almost constant) constitutes the clear 
signature by which the future experiments may show unambiguously that there is some 
non stationary process in the Sun and, consequently, in other stars as well, provided that 
uncertainties of the predicted values are less than the expected changes of the flux. By 
the present time the main contribution to these uncertainties gives the one for the cross 
section of the thermonuclear reactions which is around 10% [9J. The parameters used 
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by the standard solar model (the solar temperature et al) have a lesser impact. It has 
been shown in [10] that there is a correlation of F 13 and F 15 with F 8 so that the ratio of 
F13 (-^15) to Fg has less uncertainty than the fluxes themselves. The total flux of boron 
neutrinos has been measured precisely in the SNO experiment [11]; this makes possible 
to reduce the contribution of this uncertainty. Another important issue is the intensity 
of the mixing, i.e. the exact value of k. In the paper [6] the results were obtained for the 
case with a certain value of a mass transport coefficient (k = 10~ 10 ?/r _1 ) and for a certain 
duration T m i X of the mixing on the final stage of the solar evolution. We were tracing 
how the various parameters of the solar model have been changed with the duration 
T mix and compared the found changes with the limits established by helioseismology. It 
was shown that mixing with k = 10 _10 yr _1 and the duration T mix = 5 • 10 8 yr is within 
these limits and this can be used as a partial solution to find the general one. It means 
that a mass 0.05 of the internal zone of the Sun can be transported to and from the 
central region of the Sun during the final phase of the solar evolution of the duration 
T m ix = 5 • 10 8 yr. This indicates the intensity of the mixing process which does not 
contradict to present data of helioseismology. 

2. Detailed investigation of the case of mixing 

Here we investigate this case in more details varying freely k and T mix . On Fig. 1 
(a-d) the contours are shown for the standard relative uncertainties 1, 2 and 3 S (which 
correspond to the uncertainty 0.5%, 1.0% and 1.5% for the mean molecular weight in 
the core of the Sun) and for the abundance of hydrogen, helium, carbon and nitrogen on 
the surface of the Sun. Here 5 is a relative uncertainty which corresponds to a standard 
deviation a. 

One can see from this figure that the most stringent limit comes from the mean 
molecular weight /i in the core of the Sun and is only mildly affected by the abundance 
of hydrogen, helium and nitrogen on the surface of the Sun. This result is a good 
illustration of the power of helioseismology: what concerns the hypothetical scenario the 
structure of the center of the Sun turns out to be more crucial than what is observed 
on the surface of the Sun. Figure 2 shows the plots for the increase of the flux of F 13 
(a) and decrease of F 15 (b) relative to the ones given by a standard solar model. 

One can see that substantial increase of -F13 can be gained by mixing while -F15 does 
not experience any noticeable change. The mixing can change also the distribution of 
3 He across the radius of the Sun and, consequently, the flux of 7 Be neutrinos. The new 
set of differential equations analogous to ([1]) can be written to find these changes. 

' dX^Hj/dt = ~X 2 { 1 H)l{a 11 + a' u ) + X 2 ( 3 He)±a 33 - 

-X( 3 He)X( 4 He)±a 3 4 + HXq^H) - X^H)) 
dX( 3 He)/dt = -X 2 ( 3 He)\a 33 - X( 3 He)X( 4 He)\a 34 + 

+X 2 ( 1 H)l(a 11 +a' 11 ) + k(X ( 3 He)-X( 3 He)) () 
dX( 4 He)/dt = X 2 ( 3 He)la 33 + X( 3 He)X{ i He)\a 3A + 

+k(X Q (*He)-X( i Hej) 
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Figure 2. The plots for /i in the core of the Sun and the ones for the increase of -F13 
(a) and decrease of -F15 (b). 

where = pN& < AiAj > is the (Aj + Aj) reaction nuclear rate [3, [8]. The flux F7 of 
7 Be neutrinos on the Earth can be written as 

p3 AT 1 /■ 

®nuCBe) = Q 3 [ A A j2 J r 2 p(r)X 3He (r, t)X iHe (r, t)a u (r)dr (3) 
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where R Q - solar radius, L - Sun-Earth distance. 

One can see from Fig. 3 that the mixing compatible with the results of 
helioseismology can change the flux of 7 i?e-neutrinos by less than 2%. 
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Figure 3. The plots for [i in the core of the Sun and the ones for the increase of F-j. 



3. Conclusion 

One can see from this analysis that a small freedom is left now by helioseismology 
for a mixing of a solar matter. For periods comparable with the age of the Sun the 
mixing can be realized only with k < 10~ 10 j/r~ 1 . For greater values of k the mixing can 
exist only as the relatively short periods in the Sun's history. However, the remarkable 
signature of the mixing is the increase of the flux F 13 while all other neutrino fluxes are 
practically left unchanged. This conclusion is of general character in a sense that it does 
not depend upon the specific kind of mixing in coordinates of the internal structure of 
the Sun. From the experimental point of view it is very interesting to measure precisely 
the fluxes F 13 and F 15 . It has a discovering potential on the non stationary process in 
the Sun and other stars as well. 
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